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Abstract 

\ We study by means of renormalization group techniques the effect 

\Q ' that on the two-dimensional electron liquid may have the van Hove 

ON . singularities observed experimentally in the copper-oxide supercon- 

, ductors. We find significant deviations from Fermi liquid behavior, 

that lead to the appearance of an unstable fixed point in the renor- 

i ■ malization group flow of the effective coupling constant. Besides the 

attenuation of electron quasiparticles already known on phenomeno- 

q | logical grounds, our approach is able to explain the reduction in the 

O ■ dispersion of the band as well as the pinning of the Fermi level near 

t> . the singularity, as observed in the photoemission experiments. 
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1 Introduction 



During recent years there has been great interest in finding theoretical mod- 
els displaying the so-called "non- Fermi" liquid behavior 0]. This has not been 
a trivial search, mainly due to the fact that the Fermi liquid framework is 
quite robust in the context of interacting electron systems. The Fermi liq- 
uid concept based on Landau's theory acts as the most valuable paradigm 
in the description of real three-dimensional metals @. The relevance of such 
paradigm is not constrained to situations in which the interaction strength 
is weak and goes beyond the range of applicability of perturbation theory. 
For some reasons, not quite well understood until recently, there are not 
many ways in which an electron system, in spatial dimension D > 1, may be 
destabilized out of the Fermi liquid picture. Some of the most outstanding 
real examples in which that happens are furnished by the copper-oxide com- 
pounds that become superconductors at high transition temperatures. In the 
metallic regime of these materials, some of the most prominent experimen- 
tal properties (low-temperature resistivity, specific heat) have an abnormal 
behavior that cannot be explained by Fermi liquid theory ||. Thus, for the 
copper-oxide compounds it is not only a question of explaining the unusual 
superconducting transition temperatures, but also of understanding the un- 
conventional properties of the normal state — what may serve furthermore 
to uncover the superconducting mechanism. 

The most modern approach to a general theory of interacting fermion sys- 
tems has put forward the idea that the Fermi liquid describes a universality 
class in two and three spatial dimensions^, [5|], much in the same spirit that 
the Luttinger liquid is the class for a wide set of one-dimensional electron 
systems 0. In order to apply the renormalization group methods to many- 
body fermion systems, the main achievement has been a careful analysis of 
the interactions on the grounds of the usual classification of the perturbations 
of a statistical theory 0, |5[] — that is, relevant, irrelevant and marginal oper- 
ators. Unlike for the one-dimensional fermion systems, this is not straight- 
forward in higher dimensions, since the interaction involves in general exci- 
tations attached to different points of the Fermi surface. Thus, one has to 
talk properly of the existence of a manifold of coupling constants, not just of 
a finite number of them. Under very general assumptions, like the regularity 
of the interaction potential in that manifold of couplings, it turns out that 
there are essentially two different marginally relevant perturbations of the 
fermion system, corresponding to the superconductivity instability, which is 
always present, and to the charge density or spin density instability, which 
arises when the Fermi surface has the so-called nesting property. The rest 
of marginal interactions translates into the familiar F function of Landau's 
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Fermi liquid theory[|J. Then, in this setup there is little room left to the 
description of additional deviations from Fermi liquid behavior, which may 
be conceived only at the expense of relaxing some of the initial assumptions. 

Confronting such theoretical framework, there is clear experimental evi- 
dence that most of the copper-oxide superconductors have normal state prop- 
erties incompatible with those of a Fermi liquid. The most notorious of them 
is the linear resistivity as a function of the temperature (at the optimal doping 
for the superconductivity), while a quadratic dependence is the only behavior 
that can be accounted for in Fermi liquid theory 0. On the other hand, the 
origin of the superconductivity in the cuprate oxides remains unexplained, 
although there is strong experimental evidence in favor of a mechanism based 
on electron-electron interactions. The electrons responsible for the conduc- 
tivity in the normal state, and for the appearance of the superconducting 
condensate, are supposed to be localized on weakly coupled planes. The exis- 
tence of important correlation effects has been used to propose, in the context 
of prototypes like the Hubbard or t-J models, a novel two-dimensional elec- 
tron liquid state, that should be responsible both for the unusual normal state 
properties as well as for the high-T c superconductivity f|]. In the doped ma- 
terials, however, there is nowadays clear evidence of the existence of a Fermi 
surface typical of a metallic band|| [HJ. Furthermore, the experiments based 
on angle-resolved photoemission spectroscopy have shown interesting features 
of the dispersion relation, like almost dispersionless bands and extended van 
Hove singularities (saddle points) near the Fermi level ||10|, 11]]. Based on this 
experimental evidence, supplemented by phenomenological arguments, it has 
been proposed that the unusual properties of electrons near such singular 
points are the cause of the superconductivity in these systems^]. Although 
a variety of theoretical approaches have been attempted fl3|, [14|, [15 



full treatment of interacting electrons near a saddle point in a 2D dispersion 
relation remains to be done. 

The aim of this paper is to study the behavior of the two-dimensional elec- 
tron liquid near a van Hove singularity by means of a renormalization group 
approach. Due to the presence of the saddle point in the dispersion relation, 
the usual description of the fermion system breaks down from the start. In 
the standard many-body treatment as well as in the modern renormalization 
group approach, the fact that one can linearize the dispersion relation at 
each point of the Fermi surface plays a crucial role. When that is possible, 
one may disregard the particular features of the Fermi surface. Thus, in the 
usual many-body approach one does not care about its fate under quantum 
corrections, since all its points are regarded as essentially equivalent. When 
there is a van Hove singularity at or very close to the Fermi surface that is 
no longer true. Physically this can be easily understood since at the point of 
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the singularity the density of states diverges logarithmically. From a techni- 
cal point of view, the gradient of the dispersion relation at the saddle point 
vanishes, what makes unclear that the traditional many-body approach can 
be followed. Quite remarkably, however, the logarithmic divergences associ- 
ated to a two-dimensional van Hove singularity make the problem specially 
suitable to renormalization group schemes, along the lines discussed in ||, [J. 

In the following, we implement such a RG treatment of interacting elec- 
trons near a 2D van Hove singularity. The model to be studied is presented 
in the next section. Then, the renormalization of the coupling constant is 
worked out in section 3. In section 4 we discuss how the renormalized poten- 
tial influences the quasiparticle properties. The renormalization group flow 
of the model is presented in section 5. A discussion of the physical conse- 
quences is given in section 6. The main conclusions are drawn in section 
7. 

2 Many-body theory of the 2D van Hove sin- 
gularity 

A van Hove singularity is a saddle point in the dispersion relation of the elec- 
tron states e(k). In its vecinity, the density of states diverges logarithmically, 
in two dimensions, and shows cusps in three dimensions. The logarithmic 
divergence leads to a singular screening of the interactions, in the same way 
as for the ID Luttinger liquid [|T7j]. 

The copper-oxygen planes where the electrons near the Fermi level reside 
have an almost perfect tetragonal symmetry. Thus, the crystal wavevector k 
is defined in a square Brillouin Zone, such that — ir < k x ,k y < n. Experiments 
show that the Fermi surface lies close to the high symmetry points (0, ±7r) 
and (±7r, 0). Thus, we expect two inequivalent van Hove "points" at the 
edges of the Brillouin Zone, as depicted in Fig. 1 (note that the points 
previously defined, are, by symmetry, extrema of the dispersion relation). 

We can obtain a reasonable description of the low-energy properties of the 
system by expanding the dispersion relation around these points. Denoting 
them by A and B, and shifting the origin of momenta to each of these points, 
we can write, in general, 

E A> B(k)^(t±At)kl-(tT^t)k 2 y (1) 

where At is proportional to the deviation of the angle between the two sep- 
aratrices from 7r/2. 

Thus, the problem can be mapped onto a model with electrons with 
two flavors, which denote each of the two inequivalent singularities. The 
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full complexity of the interacting system is already present if we neglect 
this flavor degeneracy (which can be achieved, for instance, by means of an 
orthorhombic distortion which splits the two points). In the following, we will 
consider a system with electrons near a single singularity. Intersingularity 
(Umklapp) effects will be addressed in future work. 

For a single singularity, a suitable scaling of lengths leads to the simpler 
dispersion relation 



e(k) = t(kl - kl) (2) 

We describe the situation in which the Fermi energy is right at the level of 
the van Hove singularity. The case in which the Fermi level deviates slightly 
from the position of the singularity can be treated by means of a small shift 
in the chemical potential of the system, which just introduces a small infrared 
cutoff in the renormalization group approach. There remains the question, 
though, about the physical significance of a situation in which the filling of 
the Fermi sea stops at the level of the singularity. According to our point of 
view, the shape of the Fermi line may suffer changes due to the interaction 
in the system and it is not clear, in principle, that the Fermi level may find 
a stable location near the singularity. We will see later on that this tunning 
of the Fermi level to the singularity is actually quite natural, from the point 
of view of the renormalization group approach. 

We take then the chemical potential at \i = for the dispersion relation 
(H), so that the levels with e(k) < are filled and the levels with e(k) > are 
all empty. The Fermi line is given by the two straight asymptotes k x = ±k y . 
Of course, ir order to accomplish a proper description of the statistical system 
some bandwith cutoff has to be introduced that regularizes the formally 
infinite Fermi sea of (|^) . We implement in what follows a cutoff in momentum 
space {kxi ky) that preserves the SO(l,l) symmetry of the dispersion relation 
(0), by picking the states that lie in the region with — A < e(k) < A. Thus, A 
has the character of an energy (or frequency) cutoff. In fact, this possibility 
of controlling the energy scale of the processes in an invariant way is what 
allows us to carry out the renormalization program that we apply below. The 
behavior of the interactions as we reduce the cutoff A gives us information 
about the form of the low-energy effective theory for the system. In this 
respect, the procedure is similar to that adopted in the usual description of 
the Luttinger liquid [17] or in the renormalization group approach to Fermi 
liquid theory ||, 

The first step in this approach to the low-energy effective theory of the 
van Hove singularity is to write down a classical interacting theory that is 
scale invariant. We take as starting point the system of electrons with the 
dispersion relation (E3) and the typical four-fermion interaction of density- 



4 



times-density type. The action of the model reads, for the simplest case of a 
local interaction, 



S = fdt'd 2 r^^t(^m^^,t')+t^Ur,t'M-d 2 y )^ a (r^ 

J dt'd 2 r $>+(r,t')* CT M') $>+,(r,t')^M') (3) 



U 

"2 



ty a (r,t) is the electron field made out of all the modes that lie in the region 
determined by the cutoff, a being the spin index o =t, I, 

# ff (r, t) = j s duo d 2 k e - M *- T a a (k, u) S = {k | -A < e(k) < A} 

(4) 

The scaling of the different terms in (§) is best analyzed by passing to 
the frequency-momentum (uj, k) representation, in which the action reads 

S = / dujd 2 k (uj a^(k, uj)a cr (k, u) — t(k 2 x — k 2 ) (k, oj)a a (k, u)) 

J a 

~ J dud 2 k (p r (k,u)+ Pl (k,u)) V{k) (p T (-k,-^)+ Pi (-k,-^)) 
where p CT (k, uj) are the Fourier components of the density operator ^(r, t)^ a (r, t) 
Pa(K uj) = J duj p d 2 p a+(p -k,u p - u)a a (p, uj p ) (6) 

The case of a purely local interaction corresponds to V(k) = 1. In what fol- 
lows we will consider however the more realistic case in which the interaction 
has some finite range in real space. We will then suppose that the potential 
V(k) is constant but has a finite support in k space, so that the discussion 
can be carried out for a spin-independent interaction. 

We determine the scaling dimension of the operator a a (k, to) in the usual 
way, by demanding the scale invariance of the free term in the action. Given 
the nonlinear character of the dispersion relation, it turns out that under a 
change of scale of the momenta 

k -> sk (7) 

the frequency has to scale in the form 

uj -> s 2 uj (8) 

This different scaling of the frequency and the momentum is something nat- 
ural in the model, since to keep up with the approach to the Fermi level 
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e(k) = the energy has to be reduced as the square of the reduction in the 
momentum scale. From (|7|) and (§) we determine the scaling of the electron 
operator a a (k, u) that renders invariant the first piece of @ 

a a (k,uj) -> s~ 3 a a (k,uj) (9) 

A most important point is now to check whether by application of the 
scaling transformation ([?]), (§) and (0) the interaction strength U in the last 
term of (|5|) goes to zero, diverges or remains invariant in the infrared limit 
s — > 0. The first case would correspond to an irrelevant interaction, that 
should be neglected in the effective field theory at low energies, while in the 
second case the last term in (^j) would play the role of a relevant perturbation 
of the noninteracting theory, spoiling the renormalization group approach. 
We see, however, that the actual case is that the interaction remains marginal 
at the classical level, i.e. the coupling strength U does not scale as s — > 0. 
Under these conditions, the effect of the quantum corrections in the model 
turn out to be crucial since, though small they may be, they are able to 
destabilize the naive scaling in one or the other direction. 

The whole issue of the renormalization of the U parameter becomes of 
the utmost interest since, as follows from Shankar's renormalization group 
approach to fermion systems, the interaction in Fermi liquid theory (con- 
strained to the forward scattering channel) is not renormalizedQ. It has 
to be pointed out, however, that the case of a dispersion relation like (^) 
becomes special, in the sense that the interaction is marginal at the classical 
level irrespective of the kinematics of the scattering process, as we have seen 
above. This enables us to pose the problem in terms of the renormalization 
of just one coupling constant U, which affects evenly to all the scattering 
channels — apart from the BCS channel, which does require a particular 
kinematics. 

In what follows we adopt a statistical field theory description of the model, 
sticking to a path integral formulation based on the action ([|). The funda- 
mental object in this approach is the partition function 

z = J v 5( *'* +) (io) 

where, as stated above, we restrict the functional integration to the fermion 
modes within energies —A and A measured with respect to the level of the 
singularity. The Fermi sea that we are considering corresponds to a noncom- 
pact region in momentum space, though the number of states in the region 
becomes infinite only in the infinite volume limit. We will see, for instance, 
that the four-fermion interaction and the self-energy correction of the model 
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can be regularized by the use of just the bandwith cutoff A. The depen- 
dence of the quantum corrections on A gives us the information about how 
the different parameters of the model are renormalized when the cutoff is 
reduced from A to A — dA, what in practice is equivalent to perform the 
partial integration of a slice of high-energy modes in the path integral fllQ|) . 
The recursive implementation of this renormalization group transformation, 
studied by the flow as A — > of the scale dependent parameters, leads to the 
effective field theory governing the low-energy physics of the model. 



3 Coupling constant renormalization 

We study first the renormalization of the interaction potential V(k), which 
gets the quantum corrections shown by the diagram in Fig. 2, to second 
order in perturbation theory. The figure represents the contribution to the 
corresponding vertex function iT(k.,u), i.e. the part of the interaction with- 
out the external legs. f\ In terms of the fermion propagator (k, u) for 
each respective spin orientation, the vertex function to the given perturbative 
order iT^(k, u) reads 

* T ~(2^ ^^/AEGffq + ^S + ^Gftq.s) (11) 
The fermion propagator to be used in our model is 

^ 0) (q,^) = — 7? 2 (12) 

which takes into account that the dispersion relation cannot be linearized at 
the singularity. The most important point in the computation of flTT| ) is that 
the virtual states in the loop have to be kept within the band determined 
by the cutoff, i.e. only states with energies between —A and A are allowed. 
This bandwith cutoff has the virtue of preserving the SO(l,l) symmetry of 
the model, that is nothing but the invariance under the continuous set of 
transformations 

(13) 





The change of variables ([13]) preserves the form of the dispersion relation 



and it is an exact symmetry of the quantum theory. The implementation 

3 We do not mix up at this point the Cooper channel, which requires a particular 
kinematics, with the generic corrections to the four-fermion interaction, which is marginal 
in the model at the classical level irrespective of the values of the momentum exchange. 
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of a cutoff that does not break the SO(l,l) invariance of the model is an 
important ingredient of the calculation since, as we will see, it enables us to 
reconstruct the complete expression of iT^fe, to) from just the dependence 
on, say, k x and to. 

The loop integral in (|ll|) is then most easily performed by making a 
change of variables that introduces the lines of constant energy. In the sector 
IZi = {(^ \ q y > \q x \), for instance, we introduce the variables r and by 

q x = r sinh </> 

q y = r cosh (14) 



The region of integration in that sector is mapped to — oo < <fi < +oo, 
< r < A. By doing first the integration in the uj q complex plane, we get a 
partial contribution iT^ to fllTD 

?T (2) _ 
U 111 — 

= —i — - / dr d<pr — — ■ (4§-) 

2n 2 Jo J a uj — tk* — 2tk x r sinh cp + le + ze sgn (to — tr 2 J 

where a = arcsh((tr 2 — tk"l) / {2tk x r)) and we are computing setting u > 0, 
k x > and k y = 0. As we have mentioned before, the complete dependence 
of the vertex function on k and u can be reconstructed, in general, by using 
the SO (1,1) symmetry of the model, as long as the vertex can only depend 
on invariants under that transformation group. 

The angular integral in ( fLq) can be easily done, while the remaining 
integral in the r variable cannot be directly obtained. It is much simpler, 
however, to extract its A dependence, which turns out to be 

<r S wi 8^T k)gA (16) 

It can be checked that there is a similar contribution from the sector q y < 
— \q x \ and that the total contribution of the region with q x — q% > equals 



that of the region q y — ql > 0. We have therefore 



« log A (17) 

The presence of this dependence of the vertex function on the cutoff is a key 
result for the model, since it implies that the bare four-fermion interaction 
is renormalized as the bandwith cutoff is reduced towards the Fermi level. 
This property establishes a clear difference with the behavior of the forward 
scattering coupling constant in Fermi liquid theory, where it is shown that it 
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remains marginal to all orders in perturbation theory. In our case the cutoff 
dependence of can be easily understood from the physical interpretation 
of the diagram in Fig. 2. The logarithm of A in (|17| ) is nothing but a measure 
of the number of particle-hole excitations in which a particle below the Fermi 
level is promoted to an empty state above it. The logarithmic divergence can 
be alternatively seen as the divergence in the count of particle-hole transitions 
when the momentum transfer accross the Fermi line is made arbitrarily small. 

The complete dependence of the vertex function on k x and to can be 
obtained starting from the knowledge of its imaginary part. This arises from 
the sum of four different contributions like that in (|T5|). We get 



Imr(2) = i* + ; %i (18) 




Given that this object can only depend on quantities invariant under the 
transformations (|T3|), we may reconstruct the full expression 




Im r« ( k , u ) = —— l" + g WI _ (19) 



Moreover, we know that T(k, u) must have an analytic continuation from the 
half-line uj > into the upper half-plane of the complex variable uj 0. This 
fact, together with the knowledge of the imaginary part (|i9|), would suffice 
to determine the real part of the vertex function. Alternatively, one may 
introduce a different implementation of the cutoff, which makes the real part 
of r^ 2 ** directly available. The second version of the computation is sketched 
in Appendix I. There we apply the bandwith cutoff to the frequency integral 
in flllD , letting the remaining integrals in q x and q y from — oo to +oo. Quite 
amazingly, these can be done analytically by passing to the complex plane, 



producing the same imaginary part as quoted in (19) and the real part 

T2 



ReF^k,.) = — — log ^^ + ^ TIM ^ 



uj - e(k) 



uj + e(k) 



(20) 



It can be easily seen that the expression (|H|) corresponds for uj > to the 
imaginary part of the analytic function built from (|20|), with definite choices 
of the branch cuts for the logarithms. This shows the agreement between 
the two different implementations of the bandwith cutoff and finishes up the 
complete determination of the vertex function to the one-loop order. 

The renormalization of the vertex function under variations of the cutoff 
A has important consequences from the physical point of view. Within the 
path integral formulation defined by (|T0|), we may interpret r^ 2 ^ as a quantum 
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correction to the four-fermion interaction term. We define therefore the four- 
fermion effective interaction by 



.,, .1 U 2 (. |e(k)|A u , 

lV eff = -lU + log ' 7 + —r log 

27r t \ \lu 2 — e[k) 2 \ e(k) 



- e(k) 



cu + e(k) 



(21) 



The fundamental observation now is that, if V e ft is to play the role of effective 
interaction with which observables are computed in the quantum theory, it 
cannot depend on a particular choice of the cutoff A. This means that, as we 
reduce for instance the value of A, the value of the bare coupling constant 
U has to be adjusted so that the effective interaction of the theory remains 
cutoff independent. We recover in this way the idea underlying the original 
version of the renormalization group, in that the hamiltonian of the theory 
moves in the space of couplings as the length scale is varied UTS] . By impossing 
the cutoff independence of V e ff 



we get from (21), to second order in perturbation theory, 



d 1 U 2 

This equation governs the flow of the coupling constant as the cutoff is low- 
ered towards the Fermi level, leading to an effective theory for the low-energy 
processes of the model. We see that, in the case of a repulsive interaction we 
are interested in (U > 0), the bare four-fermion coupling constant decreases 
as A — > 0. It would seem, therefore, that perturbation theory should become 
more reliable at low energies if the starting value of U(A) was already small. 
One has to bear in mind, however, that the perturbative expansion in the 
model is given in powers of U/t, so that one has to check if the "hopping" 
parameter t bears any dependence on A to ascertain whether the flow at 
low energies is to weak coupling or to strong coupling in terms of U/t. We 
consider the renormalization of t in the next section, in the context of the 
quantum corrections to the two-point function. 



4 Self-energy corrections 

We undertake at this point the analysis of the self-energy function of the 
model, with the aim of studying the renormalization of the bare dispersion 
relation (0). The self-energy E(k, to) is given by the irreducible part of the 
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corrections to the two-point function and it is related to the full propagator 
by the equation 

s = ^i- E < 24 > 

The first perturbative contribution to £ is shown diagrammatically in Fig. 
3(a), where the dashed line represents the interaction potential V(k,uj). We 
have started in our model with a constant interaction l^(k), and for this 
reason it is clear that the mentioned diagram cannot give rise to a dependence 
on the external frequency or momentum. Under these conditions the one- 
loop diagram may only produce a renormalization of the two-point function 
by a constant, i.e. a renormalization of the chemical potential. We regard, 
however, this no- renormalization of the rest of parameters as a mere accident, 
since any slight modification of the potential (that still may keep it short- 
ranged) already introduces some nontrivial dependence V(k,uj). We will 
recall this observation later on, when discussing the renormalization of &(k). 

The first nontrivial contributions to the self-energy in our model are given 
by the two- loop diagrams shown in Figs. 3(b) and 3(c). Since we are dealing 
with a constant interaction the dashed lines can be contracted to a point, 
making the two diagrams look the same. The only difference between dia- 
gram 3(b) and diagram 3(c) is that the first comes with a relative factor of 
—2 with respect to the second due to the propagation of the two spin ori- 
entations around the fermion loop. The net effect is taken into account by 
computing diagram 3(b) for one of the spin orientations. We implement the 
same regularization adopted before, preserving the SO(l,l) invariance of the 
model. Thus we may take advantage of the one-loop particle-hole polariz- 
ability determined in the previous section. That is, at the two-loop level we 
have 

iZ {2) (p,u p ) = J diu k d 2 kG^\p-k,iu p -Lu k )T^(k,Lu k ) (25) 

The additional loop integral is taken within the region —A < e(k) < A. By 
dimensional reasons, it is clear that for u p , e(p) << A the dominant behavior 
is given by the terms proportional to u p and s(p), as long as the cutoff is the 
only explicit scale in the theory and these quantities saturate the dimensions 
of the inverse propagator. Those are, on the other hand, the contributions 
in which we are interested in, since the term proportional to u p in the self- 
energy renormalizes the scale of the electron field operator while the term 
proportional to e(p) may give rise to a renormalization of the "hopping" 
parameter t. 

For the sake of a practical computation of S(p, oo p ), we perform the loop 
integral by making a change of variables of the type (|13|) in the respective 
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sectors of positive and negative energy. In all cases, we first carry out the 
integral over the hyperbolic angle <p, which is finite provided we keep the 
external momentum p different from zero. To obtain the term linear in co p , 
we take the derivative with respect to uj p and set u p = 0. In the region 
k y > \k x \, for instance, we get 



— ReS (2) 



1 U 2 
16vr 4 ^" 



dr - 

r 



x 



-log 

U 2 



x 2 - 1 



n=±l 



x log 



X 



nil 



((x 
1 



n(l - e 2 )) 2 + Ae 2 ) c 



x+1 



8vr 4 t 2 



dr 



~ I'd* E 



X 



n(l-e 2 )) 2 + 4e 2 



x lot 



2e 



x 



n(l 



+ 



n(l-e 2 )) 2 + 4e 2 



(26) 



x 



where we use the parameter e = p x /r as an infrared regulator (we have 
set for simplicity p y = 0). The first term arises from the real part of I^ 2 ) 
times the imaginary part of G^ -* while the second arises from the reverse 
combination. The terms containing the explicit dependence on the cutoff 
from I^ 2 ) cancel out in S^ 2 \ what is consistent with the fact that they are 
a local insertion in the one- loop diagram of Fig. 3(a). Another remarkable 
fact regarding equation (|26| ) is that the infrared regulator e cannot be sent 
to zero, if we want to maintain finite the above expression. This will have 
important consequences, as it turns out that the dependence of on the 
cutoff becomes entangled, in principle, with a dependence on the external 
momentum. 

Putting together the pieces coming from the different integration regions, 
we arrive at the form of the linear dependence of £( 2 ) on lu p 



Re£( 2 )(p 



CO 



p 4tt 4 t 2 
1 U 2 



u 2 r A , i 

dr - 



2 log 2 log 



tr 2 



+ Ci 



CO 



4vr 4 t 2 



-2 log 2 log 2 A/ 



Px~P 



1/2 



+ Cl logA (27) 



with c\ ~ —5.896. We have used again the SO(l,l) symmetry to reconstruct 
the dispersion relation from the momentum dependence of the e parameter. 
The above expression shows that there is a renormalization of the scale of the 
electron field as the cutoff is varied, though at this point this effect cannot be 
assessed in a way independent of the electron state. In other words, to the co- 

1/2 

efficient c\ of the log A term in (p7|) we should add the log 



Pi - Pi 



coming 
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from the squared logarithm contribution, leading to the result of a wavefunc- 
tion renormalization depending on the electron momentum. This is what, 
in the field theory language, is called a nonlocal divergence, since it cannot 
be removed in the framework of quantum field theory by the subtraction of 
terms with power-law dependence on the frequency or the momentum. In the 

1/2 

present situation it may also appear doubtful that the log p x — pi log A 
cutoff dependence can be removed by a convenient redefinition of the cou- 
plings in the theory. However, the main difference of the condensed matter 
framework with respect to quantum field theory is that one may envisage the 
renormalization of the interaction potential and the Fermi surface, which are, 
in general, functions of the electron momentum. This approach to the prob- 
lem of interacting fermions is not usually adopted, mainly because in Fermi 
liquid theory the interaction, taken as a manifold of couplings, is not irrele- 
vant only in a very reduced number of channels. We have seen, though, that 
this type of restriction does not apply to our model, as the four-fermion in- 
teraction remains marginal no matter what is the kinematics of the momenta 
involved on it. The fact that the first nontrivial contribution in perturba- 
tion theory introduces a leading term of the type log 2 A, instead of a simple 
logarithm, points also at the solution of the problem. The absence of a 
log A contribution to first order in perturbation theory is just a consequence 
of having fine-tunned the bare potential V(k) to unity, which now appears 
clearly unnatural since any slight (scale independent) deviation from that 
form introduces significant renormalization effects. We will follow this line 
of thought in order to apply with success the renormalization program in the 
next section. 

We complete the discussion of the self-energy corrections by considering 
the contributions proportional to the dispersion relation e(p). Their com- 
putation is carried out by differenciating now T,^ with respect to p\ or p 2 , 
and then setting them to zero except where they act as an infrared regulator. 
Details of the calculation are given in Appendix II. We finally arrive at the 
result 



Re£< a >( 




( 2 _ 2 
VP* Py>^ t 



^(P)^i^ ((2 log 2 - ?) log 2 (A/ pl-pl *") + c 2 log Aj (28) 

with C2 ~ 3.839. We can apply to this expression the same remarks made 
for the contributions linear in u p to the self-energy. It is worth noting that 
the dependence on the cutoff in ( p8| ) is not the same as that renormalizing 
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the scale of the electron field in equation (pTp. This means that, besides the 
wavefunction renormalization also present in fl28|), there is a left-over that 
leads directly to the modification of the dispersion relation. However, the 
nonlocal character of the A-dependence prevents again to understand that 
correction, as it stands, as a pure renormalization of the hopping parameter 
t. This issue will be conveniently clarified in the next section. 



5 Renormalization 

The goal of applying the renormalization group program to a statistical field 
theory of the type we are studying is to remove all the dependence on the cut- 
off from observable quantities by an appropriate redefinition of the couplings 
(in general sense) of the theory. This philosophy corresponds to the point of 
view of relating different cutoff-dependent bare theories to a unique cutoff- 
independent renormalized theory |1| ^0[. In quantum field theory there are 
theorems that establish, in general, whether a given model is renormalizable, 
so that all the divergences as A — > oo can be absorbed in the parameters 
of the theory. In the context of condensed matter systems, however, the 
criteria for interacting fermion systems are more unclear, as the renormaliza- 
tion group approach has been developed in recent years and limited mainly 
to known paradigms, like that of Fermi liquid theory, and to models with 
properties common to quantum field theories. It turns out, though, that the 
number of degrees of freedom that are available in the renormalization of a 
statistical fermion system is much larger than that in a quantum field the- 
ory. In the latter the number of marginal couplings is always finite (except 
in the case of spatial dimension equal to one) while in a fermion system with 
a nonvanishing charge density the Fermi surface (or line) and the interac- 
tion potential should be taken as dynamical variables influencing one to each 
other. In other words, there are two different functions, not just parameters, 
that can be set to enable the cutoff independence of the renormalized theory. 
The situation is in that respect reminiscent of string theory with massless 
background fields. Although the string interaction is clearly nonlocal, there 
is always a choice of the background condensates to each perturbative order 



that keeps the theory conformal invariant [21 . 

Coming back to our model, we should allow for redefinitions of the dis- 
persion relation e(k) and the interaction potential V(k) in order to get rid of 
the nonlocal cutoff dependences of the self-energy in fl27Q and (|28|) . On the 
one hand, the saddle point in e(k) has a topological character, in the sense 
that it cannot be removed by any smooth deformation of the dispersion re- 
lation. Therefore, we concentrate on the possibility of a modification of the 
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bare four-fermion interaction by the presence of the van Hove singularity. 
There are two different ways in which this can be understood, namely from a 
physical and a technical point of view. On physical grounds, the fact that the 



renormalization coefficients in (|27|) and (28) diverge as p — > means that the 
renormalization process is out of control as one approaches the singularity. 
This instability can be interpreted as a signal that the physics is changing 
drastically at the flat region near the origin, in a way that our initial bare 
theory cannot encompass. On technical grounds, the instability can be traced 
back to the significant change that the quantum corrections bring about in 
the interaction potential at very low energies. This can be appreciated in the 
one-loop polarization tensor (|20D , which corrects the classical four-fermion 
interaction. The two effective interactions in the right-hand-side of fl2Tj| ) are 
the very source of the squared logarithm cutoff-dependence in the self-energy 
function, and they give the hint that similar terms have to be included in the 
bare interaction to absorb the nonlocal quantum corrections. From the point 
of view of a naive scaling, a logarithmic correction to the bare potential does 
not change the marginal character of the interaction, though it modifies in 
general the real space potential at long distances by a 1/r 2 tail. 

There is a sensible difference between the effective interactions that ap- 
pear in (pip, namely 



Vy ~ log uj 2 - e{k) 2 + iAi 
uj — e(k 



uj 

V 2 ~ -prr log 

e(k) 



uj + e(k) 



+ iA 2 (29) 



and those that can enter in a physical potential. A retarded interaction in real 
space gives rise to a potential in (uj, k) space that has an analytic extension 
from the real line to the upper half-plane of complex frequencies. As is 
well-known, though, the effective potentials that arise from the particle-hole 
polarizability r^ 2 ^ are analytic up to uj = 0. The imaginary parts of V\ and 
V 2 are even functions of uj, instead of odd functions as it should be to combine 
with the real part of respective analytic functions. For this reason, even if 
we renormalize the bare potential to absorb the nonlocal A-dependences that 
appear in the self-energy corrections, the effect of this renormalization can 
only cancel part of the effective interaction. 

The other obvious condition that is required to any admissible modifica- 
tion of the bare potential is to have a real part that is an even function of 
uj, leading to a real interaction in real space. Thus, we arrive at the conclu- 
sion that the real parts of the effective potentials V\ and V 2 produced by the 
particle-hole polarizability are just the only dimensionless expressions that 
can be thought of. The second potential V 2 generates however too singular 
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log A contributions when inserted in the one-loop self-energy diagram, that 
cancel only for the specific imaginary part arising from (and did not show 
up therefore in £( 2 )). It seems then that we are only left with the real part 
of V\ in our purpose of building up an analytic potential that may remove 
the nonlocal A-dependences of YP^ . There are actually two possibilities, de- 
pending on whether the argument of the logarithm is scaled with e(k) 2 or 
with uo 2 . We perform therefore a renormalization of the potential in the form 



UV{k,u) 




at log 



+ a 2 log 



\+o (u(u/tf 



(30) 

and demand that a\ and a 2 are such that, by insertion of V(k,uj) in the 
one-loop self-energy diagram, the nonlocal A-dependences cancel up to the 
two-loop order. a\ and a 2 are then uniquely determined, since as we al- 
ready remarked there are independent quantum corrections to the scale of 
the electron field and to the dispersion relation e(k). 

When used at the one-loop level, the potential ( |30"1) gives to order U 2 /t 2 
a contribution to the self-energy 
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- 2.153 a 2 log A 



(31) 



If we now set a\ = —14 log 2 + 9 and a 2 = —2 log 2 + 3, we see that all 
the nonlocal cutoff-dependences are removed from the self-energy function 
to order U 2 /t 2 

E = E {1) + E (2) + . . . 

v^i^r (" 6 - 893 lo s A) - £(p)^i^r (- 1 - 452 l0 S A) + O 



This amounts to accomplish half the renormalization program in the model, 
since there still remain the local A-dependences that have to be removed by 
the usual procedure of absorbing them in the scale of the bare electron field 
and in the bare hopping parameter. 

We require the cutoff-independence of the full Green function given by 



1 

G 



1 
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ZjftA) (u p -t(A)(pl-pl 



+Z^(A) (u p - s(p)) (6.893 log A) 

+Z^{A)e(p)~^- (5.441 log A) + O ((U/t) 3 ) (33) 

since this object leads to observable quantities in the quantum theory. In 
( |33"D Z^ 2 (A) represents the scale of the bare electron field compared to that 
of the cutoff-independent electron field 

«W(A) = 4 /2 (A)^ (34) 

From 

A-^-G- 1 = (35) 
dA 

we get the differential equations, to the two-loop order, 

K JK log Z * (A) = 6 - 893 i? (36) 

d 1 U 2 

A-y- t(A) = 5.441-^-— (37) 
dA K ' 4tt 4 t v ; 

The quantities Zy(A) and t(A) play the role of effective parameters that 
reflect the behavior of the quantum theory as A — ■> and more states are 
integrated out from high-energy shells of the dispersion relation. The equa- 



tions (pop and (|37|) have to be solved in conjunction with equation (|23f) for 
[/(A), that completes the set of coupled renormalization group equations. 

At last, we have accomplished successfully our purpose of absorbing all 
the dependences on the cutoff by renormalizing the interaction and other 
parameters of the theory. In the process we have been able to interpret the 
nonlocal A-dependent quantum corrections, that obstructed the knowledge 
about the physics close to the singularity, as a renormalization of the inter- 
action potential. This may appear a little bit surprising, since it seems that 
we had gained a complete predictability about the behavior of the theory 
at arbitrarily small energies. However, this may not be necessarily the case. 
The potentials in (|30| ) that we have introduced to renormalize the interaction 
grow large at very small values of uo 2 or e(h) 2 . Even if our renormalization 
group approach predicts that in certain regime the effective coupling constant 
U(A)/t(A) becomes small, so that perturbation theory may become reliable, 
there cannot be an absolute statement about the weakness of the interaction, 
since at very small frequencies or very large distances the interaction (|3C|) 
enters necessarily the nonperturbative regime. It seems, therefore, that we 
have not obviated yet the existence of some limitation in the close approach 
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to the singularity, at least while we do not get information about the form of 
the successive renormalizations of V(k, uf) to higher orders in perturbation 
theory. 



6 Physical properties 

The solution of the coupled set of renormalization group equations fl23|), (|36|) 
and (|37|) gives, in the limit A — > 0, the behavior of the statistical field theory 
at low energies. In principle, the nontrivial renormalizations of the wavefunc- 
tion of the electron, the hopping parameter and the four-fermion coupling 
constant are quite relevant results, since none of them have a counterpart 
in Fermi liquid theory ||, ||. The flow of the scale of the electron field in 
(p6|) , for instance, is such that it reflects the suppression of the weight of 
fermion quasiparticles at the Fermi level. This agrees with the conclusion 
obtained in the context of a purely phenomenological model of the copper- 
oxide superconductors^]. The flow of the hopping parameter in (|37| ) gives 
information about how the dispersion relation is renormalized in the low- 
energy limit. From that equation we may already draw the conclusion that 
e(k) has in any event a tendency to become flatter near the singularity. This 
behavior is also in agreement with experimental observations [[TP], ITT . 



The main issue, however, is the determination of the coupling regime to 
which the model flows at low energies. This is not just given by the flow of 
the four-fermion coupling constant U, since the parameter that enters in the 
power series expansion of perturbation theory is U/t. Thus, although U(A) 
decreases as A —>■ 0, t(A) also does and the consideration of the behavior of 
U/t requires further detail. 

From (|23| ) and (^) we obtain^ 

d XJ 1 (U\ 2 5.441 /[A 3 

A 1K7 = 2^ (?) " (?) (38) 

We regard the expression at the right-hand-side of (|38| ) as the beta function 
of the effective coupling U / 1, whose zeros correspond to fixed points of the 
flow (see Fig. 4). Apart from the fixed point at the origin, another fixed 
point shows up at a value of the coupling {U /t) critical ~ 27r 2 /5.441. One may 
question whether this point is a sensible feature of the model, as it may not 
lie in the region in which perturbation theory can be trusted. We think that 
it is realistic to consider that the perturbative expansion is more precisely a 

4 It may seem that in (^8|) there is a (U/t) 3 term missing from coupling constant renor- 
malization, but this contribution belongs to the RPA sum already encoded in the term at 



the right-hand-side of (E3 
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power series in the parameter U / (2n 2 t) , taking into account the phase space 
factors. Thus, it is very likely that the nontrivial fixed point that we obtain 
within perturbation theory is the evidence for a fixed point that remains at 
the nonperturbative level. A plot of the beta function (3(U/t) in the right- 
hand-side of (|38|) is given in Fig. 4. The points with (3 > have a flow 
that leads towards the origin as A — > while those with (3 < flow towards 
increasing values of U/t, so that the nontrivial fixed point at (U/t) cr iu ca i is 
unstable in the infrared. The solution of (|3^) can be given in closed implicit 
form 

U - U " 11 " (39) 



* 1 " 5*§ '°S (&) + «t >°S (f ft) " < '°S 



1-aU/t 



1-aUo/to 



where a ~ 5.441/(27r 2 ). It seems, therefore, that the system has at least two 
different phases. One of them is within the weak coupling regime, in which 
the theory becomes asymptotically free in the low-energy limit. The other 
phase goes into the strong coupling regime, leading presumably to different 
physical properties that cannot be studied within perturbation theory. 

Next, the behavior of t(A) can be easily understood in the following way. 
From equations ( p3| ) and ( |37D one may notice that the flow in (t, U) space is 
given by straight lines 

dt 5.441 v ; 

Then, it is clear that if we start with some initial value Uo/to < 27r 2 /5.441, 
the hopping parameter t(A) will decrease to some finite value as A — > 0. On 
the other hand, if Uo/to > 27r 2 /5.441 the running hopping t(A) will have 
no limit in its approach to zero, though our perturbative renormalization 
group methods become unreliable before reaching that value. We observe 
anyhow a definite trend of the hopping parameter towards decreasing values 
in the process of renormalization, which may persist in the strong coupling 
regime. This translates, in turn, into a tendency of the dispersion relation 
to become flatter around the singularity and, strictly speaking, a value of 
t(A) ~ would mean that the dispersion relation e(k) = t(k^ — ky) turns 
into a higher degree singularity. Quite remarkably, the experimental studies 
of copper-oxide superconductors by means of angle-resolved photoemission 
spectroscopy show very flat regions of the dispersion relation near the Fermi 
level PH|, n|. Similarly, an almost dispersionless region arises from numerical 



studies of the dynamics of a few holes in the half-filled two-dimensional square 
lattice |22[]. All this evidence points in the direction that the correspondence 



with the physically relevant situation for the cuprates may take place from 
the strong coupling phase of our model, where the renormalization of the 
dispersion towards a flat band would persist at the nonperturbative level. 
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Regarding the renormalization of the electron wavefunction, we find again 
two different situations depending on whether we stay in the weak or in 
the strong coupling phase. In the weak coupling phase, we can integrate 
equation (^) in a straightforward way taking the initial data in the interval 

(0, (U ft) critical) 

Znf(A) = e 6 - 893 /( 47r4 ) J A d p/p u{ P ) 2 /t{pf ( 41 ) 

As long as the running effective coupling U(A)/t(A) goes to zero when A — > 
according to (|3^), we can only expect a very weak reduction in the scale of the 
electron field at low energies. Thus, the attenuation of quasiparticles should 
not differ much in this case from that in Fermi liquid theory. The situation 
is completely different in the strong coupling phase, since there the scaling 
dimension of the electron field is actually changed from the noninteracting 
value. The anomalous dimension is most easily computed in the particular 
instance in which we sit at the critical coupling (U/t) cr i tica i. Applying now 
( |36| ) we obtain 

Z*(A) = A^/(^))Wt)l itical (42) 

By renormalizing then the model to the Fermi level A — > 0, the attenuation of 
the electron quasiparticles becomes complete. This feature is similar to that 
proposed in Ref. || in a phenomenological approach to the copper-oxide 
superconductors. In our framework, however, it acquires a more precise 
meaning. The renormalization effects are so strong that they change at 
very low energies the form of the electron propagator. The relation between 
the unrenormalized and the renormalized (A-independent) electron Green 
function that follows from (|42|) 

G(k, u) = A( 6 - 893 /( 47r4 ))( c/ /*)critiooi G (k, u) \ A (43) 

leads to the scaling behavior 

G(pk, p 2 u) = p- 2 +( 6 - 893 /( 2 - 4 ))(^)^-^G(k, u) (44) 

This anomalous scaling is another way of looking at the breakdown of the 



quasiparticle picture. The equations fl42|) and fl44|) are only correct at the 
critical point (U/t) cr iu C ai, that is an unstable fixed point. Therefore it remains 
as an open question to what extent the breakdown of the Fermi liquid picture 
persists when leaving the critical point in the direction of the strong coupling 
phase. 

We finish our discussion of the predictions that emerge from the renor- 
malization group approach by addressing the question of the position of the 
Fermi level in the model. The results showing a clear deviation from Fermi 
liquid behavior depend on the coincidence of the Fermi level with the energy 
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level of the singularity, as it is the divergent density of states at that point 
the source of the unconventional behavior. Actually, the renormalization of 
the observables considered above still takes place to a great extent provided 
that the energy difference between the Fermi level and that of the singularity 
is much smaller than any other energy scale in the system. It arises the ques- 
tion, however, of whether this close approach of the two levels is something 
natural, or rather any slight perturbation may send the singularity away 
from the Fermi level. This problem has been studied in Ref. p3| from the 
point of view of a system with variable number of particles and fixed chem- 
ical potential. The conclusion reached there by a standard renormalization 
group analysis is that as A — > the running Fermi energy is captured by the 
singularity. This feature is at odds with the unstable behavior of the Fermi 
energy in Fermi liquid theory, where only a repulsive fixed point is found for 
finite values of the chemical potential [|J — this fact does not have however 
physical relevance, as long as in closed systems the stability of the Fermi 
level is ensured by the Luttinger theorem. 



Rather than reproducing the analysis of Ref. [23], we give here a simpler 
but more physical explanation of the pinning mechanism. The argument 
is suited for a system that may receive extra charge from a reservoir (as 
is the case of the copper-oxide compounds) and does not have therefore 
fixed number of particles. The number of particles in the reservoir is much 
larger than that in the system, so that the flow of particles into the latter 
is controlled by the chemical potential fi of the reservoir. This quantity acts 
as a kind of external "pressure" of particles. However this external pressure 
is not seen in all its magnitude in the two-dimensional system, since it is 
partially counterbalanced there by the repulsion between the electrons. The 
competition between the two effects is what sets the value of the Fermi level 
Bp in the system. By using a simple one- loop approximation we may express 
the correction to the chemical potential as the frequency and momentum- 
independent contribution of the Hartree diagram 

e F « n + i—^ { duo d 2 kG {0) (k, u) (45) 

In terms of the density of states n(uj) we may write this expression in the 
form 

Ef = H / du n{uj) (46) 

47T J 

We may obviate the influence of the particular shape of the dispersion relation 
by differenciating (|6|) with respect to /x, so that 

dfi An d[i 
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Then, we arrive at the differential equation 

dep f 



dfj. 1 + T=n(e F ) 



(48) 



which expresses how the filling level changes in the system under variations 
of the external chemical potential. Equation ( f3£| ) shows that when the Fermi 
energy e F is very close to a level with a divergent density of states, as is the 
case of the van Hove singularity, the Fermi level is very weakly influenced 
by changes in the chemical potential of the reservoir. Of course, this does 
not prove that the level of filling has necessarily to fall close to the van Hove 
singularity, but the simple argument outlined gives a notion of the stability 
of the model with the particular filling considered in this paper. On the other 
hand, with the methods employed in Ref. it is possible to show that the 
renormalization of the Fermi energy e F keeps up with the integration of the 
high-energy degrees of freedom. The result leads to the dynamical effect of 
pinning the Fermi level to the van Hove singularity as A — > 0. Quite remark- 
ably, this feature is in correspondence with the experimental observation of 
many hole-doped copper-oxide compounds, for which the Fermi level is found 
very close to a pronounced peak of the photoemission spectraflOl. 



7 Conclusions 

We have shown that a model of interacting electrons near a 2D van Hove 
singularity can be treated by means of renormalization group techniques. 
These have found recent application in the general description of interacting 
fermion systems j|, ||, which require specific treatment given the distinctive 
character of the many-body ground state, that is, the existence of an extended 
Fermi surface. One of the main accomplishments of Refs. @, H has been 
the comprehension gained of Fermi liquid theory and its perturbations. In 
fact, under very general conditions like the isotropy of the Fermi surface and 
the short-range character of the interaction, it appears that the interacting 
fermion system falls necessarily into the Fermi liquid universality class. The 
problem that we have addressed, however, does not adhere to that well- 
established picture, as the anisotropy of the dispersion relation demands 
an appropriate treatment. As a result of that, we have found a nontrivial 
renormalization group flow of the couplings in the infrared. This flow is 
incompatible with conventional Fermi liquid theory, making the model one 
of the few systems beyond one dimension where the breakdown of a Fermi 
liquid is explicitly manifested. 

We find that, for sufficiently weak coupling, the interaction strength is 
renormalized to zero, which is an attractive fixed-point. This behavior is 
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easy to understand, as it corresponds to the standard picture in which the 
repulsive interaction is screened at long distances by the particle-hole polariz- 
ability of the medium. On the other hand, our results suggest the existence of 
another unstable fixed-point, and of a strong coupling phase beyond it. This 
new phase arises as a consequence of the tendency of the dispersion relation 
to become flatter near the van Hove singularity. As long as the hopping pa- 
rameter t is renormalized towards zero in the infrared, the effective coupling 
constant U / 1 flows to larger values in the new phase. The reduction of the dis- 
persion of the band near the singularity is in remarkable correspondence with 
the recent experimental results from photoemission studies[Tt], [TT|. There- 
fore, there are solid reasons to believe that the new phase predicted in the 
renormalization group approach should be a sensible feature of the system. 
The precise knowledge of the nature of the strong coupling phase would re- 
quire, though, different methods to that employed in this paper. Preliminary 
studies of finite clusters with next-to-nearest-neighbor interactions (exclud- 
ing therefore a spin density wave instability) seem to indicate a trend towards 
ferromagnetic order [^4j, what is consistent with the development of a very 
flat band. 

The attenuation of the electron quasiparticles that we have found in our 
model is in agreement with a phenomenological approach to the copper-oxide 
superconductors 0]. We have seen that this property arises from the electron 
wavefunction renormalization, which contains in general nonlocal contribu- 
tions that can be studied rigorously within our renormalization group ap- 
proach. The logarithmic corrections to the scale of the electron field are also 
consistent with another known property, that is, the abnormal behavior of the 
quasiparticle lifetime (~ in the model of the van Hove singularity^]. 
A standard Kramers-Kronig relation applied to the expression (|32"D leads ac- 
tually to an imaginary part of the electron Green function linear in u, in 
contrast to the conventional quadratic behavior in Fermi liquid theory. 

Finally, we have seen that the chemical potential itself also shows a non- 
trivial scaling. At low energies, the Fermi level tends to be pinned near 
the van Hove singularity. This flow of the chemical potential, unusual in 
condensed matter systems, may explain the existence of a whole family of 
compounds where experiments show that the Fermi level is close to a van 
Hove singular ity[|T0||. 

We have not explored here the role played by the intersingularity inter- 
actions, which may be important for the mechanism driving to the super- 
conducting phase. In fact, for a sufficiently large coupling between the two 
singularities an attractive channel may open, leading to pair binding with d- 
wave order parameter ||23||. It is worth noting that the Cooper channel shows 
a log 2 A dependence on the cutoff, in marked difference from the log A de- 
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pendence usual in the analysis of the Fermi liquid. In the context of our 
renormalization group approach to the model, that instability appears for a 
particular kinematics in the scattering process, corresponding to the channel 
of bound state formation, and it is therefore independent of the renormal- 
ization of the coupling constant accomplished above. It would be interesting 
then to undertake the analysis of the superconducting instability along the 
lines followed in the present paper, to elucidate whether a purely electronic 
mechanism may be responsible for the superconductivity of the copper-oxide 
compounds. 
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Appendix I 



The one-loop vertex function iT^ in ([TT]) can be computed by letting the mo- 
mentum integrals from — oo to +00, and carrying them out over the complex 
plane. In the region with io q > and for uj > 0, for instance, we have 



I+(u,k x ) = / duj q dq x dq y 

Jlu„>0 J 



" 9 >o J 00 + uo q -t [(q x + k x ) 2 - q 2 j + ie sgn(w + oo q 

1 



uj q - t{ql - q 2 ) + ie sgn uj q 
~ I du q [ dq y - 



t Ju; q >o J uj — tk 2 x - 2tk xX ju) q /t + q 2 + ie + ie ^cu q /t + q 2 + ie 

duj q [ dq y 

J -uj -tkl + 2tk xy J(uj + uj q )/t + q 2 + ie + ie 

(49) 



[uj + uj q )/t + q 2 + ie 



after picking up the residue of the two poles on the upper half-plane of the 
q x complex variable. The integrals over q y can be done by deforming the 
contour of integration in each case around a pole and a branch cut on the q y 
complex plane. The sum of all the contributions gives 

m rNitklP) arcsinwl - ljx - § 

I + (u,k x ) = — -f(u,k x ) / dx- 



t ' Jl yjx — 1 
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+—g(uj,k x ) / dx 



A/(tk%g 2 ) arcsuwl - \jx — § 



t Jl y/x — 1 
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— T9{u,k x ) \ dx 



>/( tk x9 2 ) arcshwl/x — 1 



t JO \fl — X 

2 2 

-^-\f(oj,k x )\ + ^\g(oj,k x )\ (50) 

p / * \ u-tk 2 x u + tkl 

where f(u, k x ) = ^ fc2 and g(uj,k x )= — 

The part J_ of the one-loop integral for u; 9 < can be obtained in similar 
fashion, with the result that 



%7X 

I-(uj,k x ) = —f(u,k x ) J 



) arcsh-i/— \jx 
dx — 
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-™a(u> k ) f° j / rCsh V^ 
r u}/(tklg 2 ) arcshwl/x — 1 



X 



ITT r u l\™x9 ) i 

+—g(uj,k x ) / dx- 
t Jo 

2 2 

-^-\f(u,k x )\ + ^\g(u,k x )\ (51) 

Both expressions (|50|) and (|5l| ) combine to produce a well-defined exten- 
sion to the whole (k x , k y ) plane realizing the SO(l,l) invariance of the model, 
as quoted in the main text. 
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Appendix II 



In order to compute the terms of Re S^ 2 ^(p, uj p ) linear in e(p), we divide the 
momentum integral in (k x , k y ) space in the four different sectors of positive 
and negative energy separated by the Fermi line e(k) = 0. In the region 
k y > |fc x |, for instance, we differenciate with respect to p 2 and set the ex- 
ternal momentum to zero (whenever it is not needed to regulate an infrared 
divergence). Thus we obtain 



d 

dpi 



Re£g(p^0,0) = 



1 u'- 



16tt 4 t 



/•A I roc x - nil - e 2 ) + In 

\ dr - ax } — ^^^^^=^^^^^= 

J rJo nt±i ^((x-n(l -e 2 )) 2 + 4 



47T 3 



X 

- 

r 



- log \x 2 

dx ]r 



+ x log 
x — n(l 



2 )) 2 + Ae 2 Y 



x + 1 

r 2\ i 



2n 



x log 



±i ^((x-^l -e 2 )) 2 + 4e 2 ) ; 
2e 



8tt 3 
1 

47T 3 



- n(l - e 2 ) + v /(x-n(l-e 2 )) 2 + 4e 2 

! im r ( 2 ) 



J r Jo 



dx 

n=±l 



X- 



47T 3 



/A | /■oo 
dr - dx 
r Jo ^ 



e2 )/(* - 


^(l-e 2 )) 2 + 4e 2 


x — n(l 


- e 2 ) + 2n 


(x — n(l - 


-e 2 )) 2 + 4e 2 


1 




2 ) + 


-n(l-e 2 )) 2 + 4e 2 




1 



Im r {2) 



<j/|e(k)|=x 



w/|e(k)|=x 



Im r (2) 



w/|e(k)|=x 



X 



x-n(l- e 2 ) + v /(x-n(l-e 2 )) 2 + 4e 2 



im r< 2 ) 



w/|e(k)|=x 



where e = p x /r. It can be checked that, when the contributions from the four 
different sectors in (k x ,k y ) space are put together, all the 1/e 2 terms cancel 
out. The leading contribution in loge from the integrals in the x variable 
may be evaluated analytically. The contribution independent of e is more 
difficult to extract, and it has been obtained by numerical computation of 
the above integrals, yielding the result quoted in the main text. 



(52) 
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Figure Captions 



Figure 1: Contour map of the dispersion relation for a two-dimensional 
square lattice with next-to-nearest-neighbor hopping. 
Figure 2: One-loop order correction to the interaction potential. 
Figure 3: Different self-energy diagrams to the two-loop order. 



Figure 4: Plot of the beta function at the right-hand-side of (|38 
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Figure 




Figure 3: 



